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Unstructured Low-Mach Number Viscous Flow Solver 


• Navier-Stokes equations(2-D) 


Conservation law form in terms of primative variables 

p 

substitute p = — 

Ri 


• Cell centered finite volume discretization 

• Implicit delta formulation written as: A it =~E 



Navier-Stokes equations nondimensionalized 
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Unstructured Grid Generation 

Delaunay Triangulation 
• Bowyer’s Algorithm 

Grid refinement based on aspect ratio, area, 
circumcircle radius 


Connectivity 



Grid Code Output: Geometry, Connectivity 

• Node point x, y coordinates 

• Cell nodes, cell faces, face cells 
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Preconditioning(N-S eq. 1-d) 
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Boundary Conditions 

• Implicit treatment 

• Solid wall specified as viscous no-slip or inviscid tangency. 

• Symmetry and periodic boundaries are treated through 
connectivity 


Characteristics of Preconditioned System(N-S eq. 1-d) 
Find characteristics of Ap 1 A* instead of A x . 
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when letting m -> o. Preconditioning gives a finite instead of 
infinite value of x. 



Sparse Matrix Solvers 


• Point Gauss-Seidel scheme 

• Point block Gauss-Seidel scheme 

• Conjugate gradient like method(SITRSOL) 

used to solve 

A)?=b > 


Point Gauss-Seidel Scheme 

• Every element of matrix except diagonal of block 
moved to RHS 

• Prone to divergence with poor initial conditions 

• Very sensitive to lack of diagonal dominance 

Point Block Gauss-Seidel Scheme 

• All blocks except diagonal block moved to RHS 

• Uses LU decomposition to the remaining matrix 
equation 

• More robust than the point G-S scheme 


A grid coloring scheme was used to vectorize the 
Gauss-Seidel method since it suffers from recurrence. 

The four color theorem was used to remove the recurrence 
from the convective terms. Recurrence remains in the 
viscous terms but doesn’t seem to affect the convergence 
rate. The coloring scheme was done by sweeping the 
computational cells twice. 


Conjugate gradient like solver(SITRSOL) 

• Iterative solver based loosely on the conjugate 
gradient method 

• Several iterative methods are available for 
solving non-symmetric positive indefinite 
sparse linear systems 

• Bi-conjugate gradient method 

• Generalized minimal residual method 

• Generalized conjugate residual method 

• The incomplete LU preconditioner was used 
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Results 


• Bump on Wall 

• Developing Channel Flow 

• Sudden Expansion 

• Periodic Tandem Circular Cylinders in Cross Flow 

• Four Port Valve 
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Conclusions 

• Grids can be generated about complex geometries 

• Diagonal block Gauss-Seidel solver more robust than 
point diagonal Gauss-Seidel version of solver 

• Coloring scheme allowed the vectorization of the 
implicit Gauss-Seidel solver 

• Sparse iterative solver(SITRSOL) allowed a much 

larger time step than Gauss-Seidel(ran 2 to 2.5 
times faster) — • . — 

• Temporal preconditioning allowed the compressible 
code to run at very low Mach numbers 
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